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Abstract. We introduce the Transchemistry Model as a dynamical exten- 
sion of the Algebraic Coalescence Rehadronization (ALCOR) model in describ- 
ing the hadronization of quark matter which is expected to be produced in 
' relativistic heavy ion collisions. Results are presented for CERN SPS NA49 

Pb + Pb experiment calculating hadron multiplicities and momentum spectra. 
The freeze-out properties of different hadrons are characterized by similar tem- 
perature, density and flow profiles in cylindrically symmetric geometry. 
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^ '. 1 Introduction 



The evolution of hadronic fireball was successfully described by the hadrochemical 
model 1^ in medium energy heavy ion collisions. This model is based on elementary 
Oh, processes in which hadron-hadron collisions lead to the production of other type 

of hadrons. The time-evolution of the hadronic system was described by a set of 
coupled differential equations in which the time-dependent gain and loss terms were 
determined by means of microscopical cross-sections for the elementary processes. 
The microscopical cross sections could be calculated in various approximations (e.g. 
^ ' considering free cross section, or medium effects or mean field, etc.), but the frame- 

5^ I work of the hadrochemistry calculation remains unaffected. 

On the basis of the success of the hadrochemistry the quarkochemistry model 
[|j was developed to describe the evolution of quark matter produced at ultra- 
relativistic energies. Here the constituents of the quark-gluon plasma became the 
ingredients of the microscopical processes and the cross sections could be obtained 
in various ways (e.g. from perturbative QCD). 

In the present paper we introduce a model in which the elementary processes 
contain the constituents of the semi-deconfined quark matter as incoming par- 
ticles and the constituents of the hadron gas as outgoing particles. In this way the 
elementary processes of the phase transition from quark phase to hadron phase can 
be modeled. This model, which is essentially an extension of the ALCOR model 
[|j may be called transchemistry, i.e. the chemical-like description of the phase 
transition. 
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The basic assumption is the presence of the semi-deconfined quark matter. This 
means that gluons aheady disappeared from the system (they decayed or were ab- 
sorbed) , the main degrees of freedoms are quarks and antiquarks in the deconfined 
phase, they are dressed and ready for hadronization. One can use various approxi- 
mations in the calculation of the cross sections for the elementary processes. Even 
more, one can include new terms into the set of the coupled differential equations 
to describe different collective properties of the system, e.g. temperature, flow or 
the strength of condensations, etc. 

If we make a linearized approximation to the above mentioned set of differential 
equations — assuming short time durations only — then we arrive to the Algebraic 
Coalescence Rehadronization (ALCOR) model as presented in Ref. |3|. The AL- 
COR model is very successful to describe the hadron yields produced in different 
high energy heavy ion collisions Q . Both transchemistry and ALCOR model yield 
particle numbers in the reaction volume V which can be compared with the experi- 
mental results. However a more complete comparison can be developed if we know 
the full momentum distribution of the out-coming hadrons. Since cr, the production 
(hadronization) cross section has its maximum when the relative momentum of the 
participants is small (see later eq.(|2^), then we may assume that the spectrum 
of the produced hadrons can be characterized similarly as that of the 
quarks. Although a direct calculation of the momentum spectra of the out-coming 
hadrons is in progress, which is based on the microscopic coalescence applied in AL- 
COR j^], in the present paper we use the approximation mentioned before. Since 
the hadrons will suffer inelastic and elastic collisions during the time-evolution of 
the hadronic fireball, one can assume that all hadrons can be characterized with 
the same local temperature (or with the same temperature profile in case of an 
extended particle system) . Because of the presence of a quark matter phase before 
hadronization one can assume that the collective flow profile will also be the same 
for the different hadron species. Furthermore, we will assume, that all hadrons can 
be characterized by the same density profile. The collective behaviour combined 
with the particle numbers obtained from transchemistry or ALCOR yield the nec- 
essary constrains to obtain the momentum spectra for all hadrons if we determined 
one of them (e.g. the spectrum of pions) from the experimental results. 

In Section 2 the basic equations of transchemistry will be summarized and we 
show, how the former version of the ALCOR model can be derived from the equa- 
tions of transchemistry. In Section 3 we display the ALCOR model, emphasizing 
its main properties. In Section 4 the hadronization rate will be introduced in the 
presence of flow which is a general result for particle systems characterized by a 
flow pattern. In Section 5 we display the particle numbers obtained from ALCOR 
for the Pb+Pb collision at SPS energy. We also calculate the fully equilibrated 
hadron gas scenario and we discuss the main differences between these models. In 
Section 6 we generate momentum distributions for pion and proton spectra and 
compare with the experimental results in the Pb+Pb collision. 
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2 Transchemistry and ALCOR 



In this section we will summarize the basic equations of transchemistry. For sim- 
plicity let us assume, that the net baryon charge of the system is zero and only 
TT mesons are produced. Thus the rate-equations for pion production from semi- 
deconfined quark matter are the following: 

I [V{tM)] = D^^\a^^v)^t) ■ n,{t) ■ n^V{t) - n^{t) ■ V{t) ■ (1) 
The rate equations can be expressed in terms of particle numbers, too. 

j^N^it) = ^F(^f^ • ^S) ■ m) - NAt) ■ r., (2) 

One obtains similar equations for the quark and antiquark numbers: 

j^N.it) = A_2^.^^(t).Ar-(t) + 7V^(<).r. (3) 



j^N^it) ^^^^■N,{t)-N^{t) + NAt)-r^ (4) 

Here ni{t) denotes the density of particle i, V{t) is the reaction volume within 
which constant density was assumed and Ni{t) is the number of particle i. The 

denotes the loss term for pions related to their 'decay' (or re-entering the semi- 
deconfined quark matter). The {'^qq'>^)(t) denotes the momentum averaged elemen- 
tary hadronization rate for one hadronic degree of freedom and D'^'^^ is the spin 
degeneracy factor for the final hadron. The temperature, density and flow profile 
for the quark distributions inside the volume V can be chosen in various ways. 
We will assume that they are identical with the hadronic ones obtained from the 
experiments. The hadronization rate can be calculated in many microscopic ways 
assuming different models: quark coalescence Nambu-Jona-Lasinio model 
or other effective quark models, e.g. Q. 

From eq.(2)-(3)-(4) one obtains that N^^ + Nq = const, and N„ + Nq ~ const., 
which is the number conservation for the hadronization process q + q <=^ tt. If the 
microscopical hadronization mechanism is different, e.g. q + q tt + tt, then the 
above equations will contain different coefficients. In its recent form transchemistry 
equations conserve not only baryon numbers, but the numbers of the constituent 
quarks and antiquarks. 

In general case the particle number production for meson M can be described 
as follows: 

-NMiqr,qj;t) = NqMN-q^{t) - NM{q^,q,■,t)^M (5) 
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In the loss term the J2 means all of the possible decay channels. We plan to apply 
the transchemistry for fast hadronization processes, in which case the deconfinement 
channels have minor role. Thus we will neglect the decay channels in the foUowings. 

If the hadronization starts at time to and ends at time to + t, then one can 
obtain the produced hadron numbers from the following integral: 

NKM,qfM + T) = ^ ■ N-,^{t) dt (6) 

In transchemistry various hadrons can be formed from the deconfined quark 
matter. Meson formation is relatively simple, because one step is enough to confine 
a quark-antiquark pair into meson-like hadron. Baryon (and antibaryon) produc- 
tion might happen in three-body interactions or sequentially via two-body interac- 
tions. In the latter case diquarks (anti-diquarks) will be produced and the time- 
evolution of their density should be followed also. The microscopical description 
of baryon production is not a well-developed field, so we use the simplest two-step 
processes to form baryons. In this case we need to follow the production of diquarks 
and anti-diquarks also, as they are carrying explicitly certain degrees of freedom in 
transchemistry. 

In case of fast hadronization (e.g. an overcooled quark matter is converted 
into hadrons [^) the volume V and the averaged rate {av) will not change too 
quickly, thus one can apply the theorem of mean value of the integral calculus. 
For one integral this transcription is exact. The first term in eq. (|^) — which 
contains the rate and the volume — remains approximately unchanged, but the 
time-dependent function of quark numbers will be substituted by their mean value 
at time t* {t^ < t* < to + r, where r << to): 

NM{q^.qfM + r) - ^ iV,. (t* ) • iV^^ (t* ) (7) 

One can approximate the particle numbers at t* by means of the initial particle 
numbers: Nq.{t*) « bq.Nq.{to) and Nq,{t*) « bq.Nq.{to). If the parameters bq. and 
bq^ were time-independent then the above approximation can be useful to simplify 
the time-dependent transchemistry equations: 



NMiq^,q,;to + r) = ^ bq^Nq^{to) ■ bq^NqAto) (8) 



^(M)(^_„).. 

V 

The complete set of time dependent differential equations of the type of eq. (||) 
preserve all particle numbers which should be preserved — as we mentioned earlier. 
However, after the above linear approximation the conserved particle numbers will 
not be preserved automatically. Therefore we have to impose the conservation 
law of particle numbers (namely quark and antiquark numbers) while keeping the 
essence of the above set of equations. In this way one arrives to the basic equations 
of the phenomenological ALCOR model, as a linearized, averaged special case of 
the transchemistry. 
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3 The ALCOR model 



The starting point of the ALCOR model [QQ is a semi-deconfined state in which 
only those quarks and antiquarks are present in local thermal equilibrium which 
form the final hadrons. We also assume that only quarks and antiquarks are present, 
so their number already accounts for those gluons which were fragmented earlier. 
This assumption is easily acceptable for the heavy quarks for which any meeting 
with their antiparticle is improbable due to their small specific density. In the case 
of light constituent quarks (counted after gluon fragmentation) the approximate 
entropy conservation during the rehadronization process leads to the conservation 
of their numbers. This picture was first formulated in the framework of the alge- 
braic recombination model ||^ . In this model the number of a given type of hadron 
produced is proportional to the product of the numbers of their constituting quarks. 
This picture was developed further by introducing gluon fragmentation We 
shall denote the number of produced u, d, s quark pairs — just before the hadroniza- 
tion — by N^^pair, Nd^pair and Ng^pair, respectively. The strangeness production 
factor, gs is defined as gs = N pair / [Nu, pair + Nd,pair)- 

The mesonic coalescence factors of the ALCOR model can be introduced by 
means of eq. (^), namely 

CM{q^,q,) = (9) 

Thus the production of meson M can be described as follows: 

NM{q^.q,) = D^^'^CMiq^,q,)■b,^N,^■b^^Ng^ (10) 

The baryonic (antibaryonic) coalescence factors can be defined and calculated as 
a simple sequential process, when baryons containing three constituent quarks will 
be produced from diquarks and quarks. The coalescence factors are the following: 

CB{a,b,c) = gB'^{CM{a,b)CM{[a + b],c) + 

+CA/(a, c)CMi[a + c],b) + Cnib, c)Cm([5 + c], a)} . (11) 

Here indices a, 6, c means the qi quarks for baryon production or antiquarks for 
antibaryon production. The factor is the baryon suppression factor. Its physical 
meaning can be connected to the sequential creation of the baryons. We assume 
that the production of a baryon from three quarks or an antibaryon from three an- 
tiquarks is a two-step process leading through an intermediate diquark formation 
in the corresponding color triplet or anti-triplet state which forms the final baryon 
later on together with a third quark or antiquark. These diquarks must, how- 
ever, be very unstable, short-lived clusters, so they may decay before forming the 
baryon. We take into account this mechanism of baryon production by introduc- 
ing the phenomenological baryon suppression parameter, gs- Mathematically one 
can introduce the factor gs similarly as the parameter hi, namely approximating 
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the time-dependent integrals for baryon production with the substitution of mean 
integral values. However, gs seems to be a well-behaving physical parameter, so 
we will determine it from the experimental results, similarly to gs- We will assume 
that gs is flavour blind. 

Now one can determine the numbers of the baryons and antibaryons: 

Nb (<Z^ , , gfc) = D^^^ Cb (<?,; , <Z, , ft) • A^,, • iVg, • hq, Ng, , (12) 

Ns{q,,qj,qk) = C^{q,,q^,q,) ■ b^^N^^ ■ b^^N^^ ■ b^^N^^ . (13) 

In the above equations the spin degeneracies are also included, D<-"^ = 25*^ + 1. 
Here we shall consider the simplest hadron multiplets conserving isospin symmetry. 
These are the spin and spin 1 meson octets together with the a and to isospin 
singlet meson states. Furthermore the spin 1/2 octet and spin 3/2 decouplet baryons 
and anti-baryons are also included. 

Since we have introduced the numbers of different mesons, baryons and an- 
tibaryons produced from the quark matter, one has to put up the final set of 
equations, which will ensure the conservation of the numbers of different flavours 
(for each quark and each antiquark flavor) : 

Nf 

H j=l 
Nf Nf 

+ E E E(l + + kk)D'^''^ CB{q^, q,,qk) ■ b,,N,^ ■ b,^N,^ ■ b,,N,, (14) 

H 3 = 1 k=j 
Nf 

= E E D^''^CM{q,. qj) ■ b^^N^^ ■ b,^N,^ 

H j=l 
Nf Nf 

+ E E E(i + + kk)D^'''>Csiq^.qrqk) ■ K^-,^ ■ • hN-,,. (i5) 

H j=lk=j 

This set of equations determines uniquely the normalization factors, bi. The 
number of independent equations in eqs.(p^,(p^) is equal to the number of inde- 
pendent 6g. and 6g. factors. Note, that these 2Nf quantities are not adjustable 
parameters of the model, but they are determined by the flavour conservation. 

Furthermore, introducing the factors bq. and bq. we have constructed a scale 
invariant description for the particle production. Namely, the physical quantities 
can be transformed in the following way: 

Ca/ {qt,qj) — » XiXjC M [qi ,qj), 

CB{qi,qj,qk) — > K Xj Xk CB{qz,qj,qk), 

% ,qj,qk} — ' ^« ^fc % (9i ,qj,qk): 

bq, > bqJXi, 

hi bqj^^ (16) 
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Substituting the results of this transformation into eqs. (|To|), (p^, (13), one 
can see that the numbers of produced hadrons, NM{cii,qj), NB{qi,qj,qk) and 
N-g(q^,qj,qf.) remain unchanged. Thus even if we miss a constant factor in the 
microscopical hadronization coefhcients, then one can obtain the particle numbers 
correctly. 



4 Hadronization Rate in the Presence of Flow 

In the presence of a flow we consider the relativistic Jiittner distribution 

/(x,p) = e-'^P-"^^) (17) 
with Ufj_{x) being the local four- velocity of the flow. Let us restrict ourselves in 



the following to a scaling longitudinal flow 1 11 . In this case the thermal averaging 
of the rate can be interpreted only locally: the earlier rate equations has to be 
reinterpreted in terms of reacting components from different coordinate rapidity 
ranges rja and r]b respectively 

\ d f 

y-^{Vnh)= driadrib{crv)abnanf,. (18) 

Clearly the total volume of the fireball increases due to a longitudinal scaling expan- 
sion like V = tttE? At/ with total coordinate rapidity extension A?7 and transverse 
radius R. Assuming a space-time rapidity plateau for a finite interval in the number 
distribution of the reacting and produced components we arrive at a modified rate 
equation 



r^Nh = -^NaNb ■ (19) 
dr ttR^ 



Here 



A* = y d{T^a - Vb) Wv) e{{pa - Pb) ■ {Xa " Xb)) (20) 

is the total rate of reactions between all possible coordinate rapidity cell pairs. The 
constraint {pa ~ Pb) ■ {xa — Xb) > reduces to the requirement that the relative 
velocity vector is oriented opposite to the relative position vector in the center of 
mass system of the colliding components 'a' and 'b'. This constraint is necessary 
because only those particles collide which have a relative velocity pointing towards 
each other. 

The inclusion of this constraint after some algebraic manipulation leads to the 
general expression of the averaged rate 

SmlmlK2{f}ma)K2{fimb) ^ ' 
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where A*^(s) — [s — {nia + rrib)^] [s — {nia — nib)'^] and G{f3^/s) is the thermal 
weight factor. In the presence of longitudinal scaling flow the latter can be written 
as an integral over relative coordinate rapidity -q = [rja — rjb) 



Ko (/3y5ch|7y|) - Ko (pV^ch\rj\ + /3^/^sh|7y|') 

G{(3V~s)^ I drjVs J_ ^. (22) 

I^VKbis) chh|sh|77| 



Investigating the integrand of eg. (|22|) one observes a finite width, Sj], in the relative 
coordinate rapidity rj. (This quantity, drj, is equal to the relative flow rapidity be- 
cause of the Bjorken scaling assumption.) The finite width depends on temperature 
T, on the particle rest masses m^, rrib, and on the considered energy scale ^/s in a 
complicated manner. 

The most uncertain ingredient of the ALCOR hadronization model is the had- 



ronization cross section in eq. (21). We will consider an analogy with the p + A 
d+ {A~l) nuclear rearrangement (pick-up) reaction leading to deuteron formation 
[ p^ . Assuming a Coulomb- like potential in the final (qq) state a fusing cross section 
can be derived from this analogy: 

2 

a = WmlV^p' " ° (23) 
(1 + [kaY) 

Here nih is the rest mass of the meson, while a = 1/ {vriabo) is the Bohr radius of the 
bound qq state in a V{r) — —a/r Coulomb potential with niab being the reduced 
mass of particles 'a' and 'b'. In our calculation we took a — 0.46 and T — 200 
MeV temperature in estimating hadronization cross sections at CERN SPS energy. 



The factor p — 0.3 fm occurring in cq.(23) accounts for the medium influencing 



the hadron formation and was taken to be equal to the Debye screening length in 
quark-gluon plasma at the above temperature. Finally k occurring in eq.(P3|) is 
the magnitude of the relative momentum vector of particles 'a' and 'b' measured 
in their center of mass system k = ■^A*j(s)/2y^. 



5 Results for Pb+Pb collision from ALCOR 

Well armed with the above theoretical considerations the redistribution of different 
flavor quarks and antiquarks into all possible hadrons can be calculated in the AL- 
COR model using practically only three parameters: i) the total number of quark- 
antiquark pairs, Ntot.pair = Nu,pair + Nd^pair + N^^pair, which can be determined 
from the measured total charged multiplicity, ii) the parameter gs controlling the 
baryon formation and iii) the strangeness production factor gs- The dependence of 
the final results on the other parameters (e.g. T, a) within their physically accept- 
able interval is small. All further results, such as the number of hyperons, kaons, 
etc. are predictions of the ALCOR model. 
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In Ref. 0] we presented our results for the S+S collisions at 200 GeV/nucleon 
bombarding energy. We took into account that the number of participant nucleons 
was measured to be -/V^^^^j^ = 51 |]l3j. We used the parameters Nf^f p^^^^, — 158.1 
with gs = 0.255 for the strangeness production, and gs = 0.04 for the baryon for- 
mation. The obtained particle numbers fitted the experimental results excellently. 

From this experience we conjectured that in S+S collision the rehadronization 
process is locally quick: the quarks and antiquarks (including the fragmented glu- 
ons) appear in hadrons according to simple kinematic rules and production branch- 
ing ratios. This finding, however, does not settle the question whether quark matter 
with any degree of collectivity has been formed or other hadronization processes 
including strings or color ropes are the source of this amount of quarks. Local 
thermal distribution of the different quark flavors and the presence of a longitudinal 
flow cannot, on the other hand, be excluded on the basis of this experimental data. 

To obtain results for A+A collisions we assume a scaling connected to the par- 
ticipant nucleon number for the produced quark-antiquark pairs: 

<ta.riV-s)=[-0^] Nglp^UVS), (24) 
Y partic 1 

where N^^^^^ are the number of participant nucleons in the A+A collision. The 
scaling exponent a may have the value a = 1, or, for more collective production 
processes one expects a > 1. Furthermore we shall use the values for gB and gs 
obtained above. 

At the extrapolation of the total number of produced quark-antiquark pairs 
from the S+S collision one has to consider an energy scaling. We will assume a 
logarithmic one: 



In our case, considering 160 GeV/nucleon collision energy for the Pb+Pb collision, 
this rescaling yields an « 4% correction for the N^^^'^.-^,. 

Now we can make prediction for the total particle production in the Pb+Pb 
collision. The number of participant nucleons is N^^^^ = 390 ± 10 obtained from 
Monte-Carlo simulations. Using the mean value and a linear scaling (a = 1) one 
obtains N^Jf^^^^ = 1164. We keep the strangeness production and baryon forma- 
tion factors fitted for S+S collision, gs = 0.255, gs — 0.04. Table 1 shows our 
prediction for the particle numbers together with the results of two thermal equi- 
librium models. Thermal I and Thermal II. In Table 2 some characteristic particle 
ratios are displayed for the three models. 
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Table 1. Particle numbers: prediction of ALCOR model for the Pb+Pb col- 
lision and the results of full equilibrium thermal models, 'Thermal I-IF. 



TDK 1 T>"U 




Inermal 1. 


Inermal 11. 


h~ 


730.41 


798.72 


839.44 




im.si 


(i()3.97 


(i()3.")4 




618.95 


620.70 


619.16 




D04.D0 


DOi .yi 


000. OL 


K+ 


84.15 


174.57 


205.68 




84.15 


177.74 


211.29 


r 


41.65 


103.10 


119.46 


K- 


41.65 


104.86 


122.69 




62.90 


140.42 


165.38 




170.90 


141.58 


129.97 


Tl 


Loo.D 1 


1 /I on 


1 QQ A A 
iOO.44 


s+ 


12.86 


24.21 


41.27 




13.63 


25.04 


42.37 




14.43 


25.89 


43.51 




68.20 


55.87 


42.37 




8.82 


10.92 


11.99 




8.89 


11.09 


12.31 




1.48 


1.56 


2.92 




25.07 


13.28 


14.57 




25.07 


12.68 


14.19 




4.18 


4.10 


8.10 




4.18 


3.96 


7.89 




4.18 


3.83 


7.68 


A^' 


20.93 


8.84 


7.89 


— t) 


5.98 


3.10 


3.91 




5.98 


3.05 


3.81 




2.20 


0.77 


1.59 
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Table 2. Particle ratios: prediction of ALCOR model for the Pb+Pb collision 
and the results of full equilibrium thermal models, 'Thermal I-IF. 



Pb+Pb 


ALCOR 


Thermal I. 


Thermal II. 




0.28 


0.23 


0.21 




0.14 


0.09 


0.11 


p-/7r- 


0.04 


0.02 


0.02 


K+/t:+ 


0.14 


0.29 


0.34 


K+/K- 


2.02 


1.66 


1.67 




0.92 


2.51 


3.90 


r/-p- 


0.83 


0.66 


0.54 




0.30 


0.16 


0.18 


S-/A" 


0.13 


0.19 


0.29 




0.28 


0.34 


0.48 




1.48 


3.62 


3.22 


p-/K- 


0.60 


0.13 


0.12 




1.64 


0.53 


0.35 



In the thermal models we assume a homogeneous system in full chemical equi- 
librium at different level. In model 'Thermal-I' the same hadron multiplets are 
assumed to be in equilibrium, which are populated in the ALCOR model, and then 
the short lived resonances are allowed to decay into the stable particles shown in 
the Tables. On the other hand, for the model 'Thermal-II' we assumed, that the 
final (stable) particles are in chemical equilibrium. In both cases the values were 
determined to get approximately the same tt"*" number as in ALCOR at a given 
temperature, T = 165 MeV. The chemical potential for strange quarks were de- 
termined to reproduce strangeness neutrality. The light quark chemical potential 
was obtained from the total baryon number of the system. 

Comparing the results from ALCOR with the results from the thermal model, 
we find that, surprisingly enough, they are not far from each other. The greatest 
deviation is about a factor of five in some of the particle ratios (see e.g. /K^ ). 

This similarity is due, probably, to the fact that, i) we have a large number of 
resonances; ii) the spin degeneracy factor appears in both models; iii) the conser- 
vation rules lead to strong restrictions. 

Besides the similarities there are systematic differences between the two models. 
The strange baryon to strange meson ratio in the thermal model is systematically 
smaller than in the ALCOR results. This can be interpreted as the effect of mass 
differences, which leads to an exp(-m/T) suppression in the thermal case, while it 
is not present in ALCOR. Thus this ratio may indicate that how strong was the 
hadrochemical evolution in the final hadronic phase. 
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6 Momentum distribution 



In the ALCOR model we have a microscopical dynamical prescription for the pro- 
cess of quark coalescence into hadrons. Furthermore we assumed a given momentum 
distribution for the constituent quarks. Thus the momentum distribution of the 
produced hadrons can be calculated. However, this is a somewhat tedious proce- 
dure. On the other hand we observe, that the capture cross section ( eq.(p3|) ) has 
its maximum at zero relative quark momentum. Therefore one may assume, that 
the hadrons, just after their production will have a momentum distribution similar 
to that of the quarks ( eq.(17) ). In this chapter, however, for the hadrons wc shall 
use Fermi and Bose distributions, and a transverse flow is added, too (see eq.(28) 
below). We assume, that, due to the elastic collisions, all the hadrons will have 
the same density, flow and temperature profiles. However, the particle composition 
will remain the same as it was at the hadronization. Thus wc have chemically not 
equilibrated system, where all hadrons will have its own chemical potential, corre- 
sponding to the partial density of this species. The value of the temperature, T, the 
transverse radius, R, the radial flow velocity at the surface and the proper time, r, 
are determined by fitting the measured rapidity and transverse mass distribution. 
However, due to our assumptions, the same set of parameters must be used for all 
type of particles. Thus the existence of such a self-consistent fit is not trivial! From 
the momentum distributions and the volume we can calculate the entropy of the 
system, too. 

In general the invariant momentum distribution will be the following: 

_ 1 f 9deg y, i^p. 

We will keep the quantum statistics, so for mesons we have K = —\, for baryons 
and antibaryons K = \. Later we will see the importance of quantum statistics 
for pions. For correct freeze-out calculation we need to introduce the Cooper - 
Frye pre-factor p'^dcTjy [ p^ . Let us consider a cylindrically symmetric geometry and 
assume freeze-out at constant proper time r. For such a freeze-out hypersurface 
one obtains 

p'^doi, — nit rdr Tcosh{y — rf) drj d(/)r- (27) 

Furthermore, let us assume a longitudinal Bjorken flow and a linear self-similar 
transverse flow, 

~ [cosh{r])cosh{Q), sinh{r])cosh{Q), sinh{Q)cos{<^r), sinh{Q)sin{(l)r)] (28) 

where the transverse flow is defined as Vt = tanh(0). The transverse radius depends 
on the proper time during expansion, namely it is R{t). 
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After trivial calculation one obtains the invariant momentum spectra: 







/ rdr 1 


drj / 


'0 Jr)^i,^ 





mtdrritdy 2Tr {hc)^ J „ 7,,„„, Jo etP""-"'")/'^ + if 

where 4> = 4>r — 4>p and p'^u^ = rrit cosh(0) cosh(?/ — r]) — pt sinh(0) cos ip. 

The calculation of particle spectra is based on the eq. (p9|), where the thermo- 
dynamical parameters depend on the radius, T{r) and fi^r). Considering Bjorken 
scaling, we neglect here longitudinal dependence, only the finite longitudinal size 
of the firebaU is included, [zmin, Zmax] = [t ' sinh?]™^, r • cosh?7„i„]. 

If we investigate the radial dependence of the particle numbers, we need to 
introduce the following distribution function. 



rdr 27r(?ic)3 7^_ '7_oo Jm Jo ^e(p-"-A^)/^ + if 

We assume that these radial distributions are similar for different particle species 
at freeze-out time t/. 

The total energy of the system can be obtained from eq.(30) by multiplying the 
integral with E — mt coshy and integrating over r. 

Before comparing the experimental results and the calculated particle spec- 
tra, let us summarize shortly the possible values of different parameters. For the 
temperature profile no longitudinal dependence was considered, the radial depen- 
dence was chosen to be quadratic or higher, T{T],r,T) = To(t) • [1 — (r/i?(r))^], 
To(r/) — 150 — 180 MeV, 7 > 2. There is a longitudinal Bjorken scaling and 
for the transverse flow we will use the parameterization vt{t,t) = tanh(0(r, r)), 
where sinh(0(r, r)) = wt{t) ■ r/R{T) and on the radial surface at freeze out 
time it is VT{R{Tf),Tf) = WT{Tf)/y/w^{Tf) + 1 < 1. The transverse radius, 
R(Tf), the minimal and maximal value of the space-time rapidity, {rjmin, Vmax) 
and the proper time of the freeze-out Tf are free parameters (zmm = t smhrjmin, 
Zmax = T sinhrimax)- Chemical potentials, /i(r) will be calculated from dN/rdr. 

Fig. 1 and 2 display the calculated n" and proton spectra at Tf = 6 fm, 
Toi^f) = 165 MeV, 7 = 2, rimin,max = ±1-75 (which means longitudinal size 
1^1 = ±16.75 /m), RriTf) = 14 /m, WT(i?(T/)) = 0.77. In Fig. 2 the full lines 
indicate Bose and Fermi distributions, the dotted line indicates results for pion with 
Boltzmann distribution. The transverse momentum distribution of pions shown in 
Fig. 2. suggests that the Bose statistics plays a role in the enhancement of the low 
mt part of the pion spectra. Such an effect was suggested earlier for the heavy ion 
collisions at Bevatron ||l6| . If we increase the freeze out volume by a factor of two, 
the pion density decreases, and the low mt enhancement disappears. 

The total entropy and energy of the hadron system were discussed in Ref. JitI , 
where it was obtained S/Nb = 30 for the specific entropy in the above parametrized 
hadron gas. 
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Fig. 1. Rapidity distributions of pions (top) and protons (bottom) from the 
extended ALCOR model for 160 AGeV Pb + Pb central collisions. The open 
squares are the measured and the reflected data points for h- particle (top) 
and '(+) - ('-^'particle (bottom) spectra ^j. The value of proton multiplicity 
at mid-rapidity is uncertain. In more recent publication [^] a small drop was 
given, but even the statistical errors are appreciable, yet. 
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Fig. 2. Transverse mass, mr, distribution of pions (top) and protons (bot- 
tom) calculated from ALCOR for 160 AGeV Pb+Pb collisions at mid-rapidity. 
The dotted line in the upper figure shows the result obtained with Boltzmann 
distribution. The open squares are the measured values for h- particle (top) 
and '(+) - ('-^'particle (bottom). 
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r(fm) 

Fig. 3. The approximately similar radial density profiles for pions and protons 
(top) and the common temperature profile, T(r) (bottom) together with the 
calculated chemical potential for protons, fiproton- The pion chemical potential 
was constant, h-k = 103 MeV. 
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We would like to point out, that the description of the particle spectrum from 
the extended ALCOR model is based on a parameterization of the flow, density 
and temperature profiles of the hadronic phase in the time interval when the hot 
hadronic matter freezes out. This parameterization is qualitatively similar to the 
ones developed in refs. |2^]-|2^, for the description of the single- particle spectra 
and the correlation function of the tt and K mesons, as has been measured by the 
NA44 collaboration. 

Both models confirm that the transverse flow profile, the local density profile and 
the transverse temperature profile |^l|-[p3| may play a crucial role when describing 
the hadronic spectra. 

Conclusions 

In this paper we introduced the concept of transchemistry for the description of 
hadronization. The ALCOR model is a good approximation to the transchemistry 
description, if the hadronization is fast. The ratios of different hadron multiplici- 
ties obtained from ALCOR and from the thermal model - based on thermal and 
chemical equilibrium of hadron resonances - are much closer to each other, than 
one could expect. The rapidity and transverse mass distribution from ALCOR was 
compared to the experimental ones from the SPS PB+Pb collisions. The agreement 
is good. 

Further we have shown that the Bose statistics may play a role in the low rrit 
enhancement in the transverse mass distribution. 

These results support the claims that in high energy heavy ion collisions a semi- 
deconfined state of the matter is formed in Pb + PB reactions at CERN SPS 
energies. 

Since it seems that in the particle composition there are only minor, but sys- 
tematic differences between ALCOR and thermal models, refined analysis tools are 
expected to play a critical role in the future data analysis: especially fits to 
the slope parameters and to the detailed invariant momentum distribution of var- 
ious kinds of detected particles may be necessary to clearly distinguish among the 
proposed models. 

Note also that the total volume and the inner dynamics characterized by the 
local temperature, fiow and chemical potential distributions are very different in 
ALCOR and in the models assuming global thermal equilibrium. The HBT radius 
parameters are known to be very sensitive to the detailed structure of the flow, 
temperature and chemical potential proflles. Thus we expect that their evaluation 
using parameters fitted to the spectrum can possibly clearly distinguish between 
ALCOR and simple thermal models. Further work is necessary to evaluate this 
parameters from the models. 
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